Competition between antiferromagnetism and superconductivity, electron-hole 
doping asymmetry and "Fermi Surface" topology in cuprates 
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We investigate tlie asymmetry between electron and hole doping in a 2D Mott insulator, and the 
resulting competition between antiferromagnetism (AF) and d-wave superconductivity (SC), using 
variational Monte Carlo for projected wave functions. We find that key features of the T = phase 
diagram, such as critical doping for SC-AF coexistence and the maximum value of the SC order 
parameter, are determined by a single parameter r] which characterises the topology of the "Fermi 
surface" at half filling defined by the bare tight-binding parameters. Our results give insight into 
why AF wins for electron doping, while SC is dominant on the hole doped side. We also suggest 
using band structure engineering to control the rj parameter for enhancing SC. 
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PACS numbers: 74.72.-h, 74.20.-z, 75.10.Jm, 71.27.-|-a 

Ever since their discovery, cuprates continue to 
pose some of the most challenging theoretical puz- 
zles [l[ in condensed matter physics. The problem is 
dominated by strong electronic correlations 
and the t J-model and its variants are believed to 
contain the essential physics of high Tc superconduc- 
tivity. In this Letter we address the following ques- 
tions: What controls the electron-hole asymmetry in 
cuprates? Why does antiferromagnetism dominate 
on the electron doped side, and superconductivity on 
the hole doped side? How can we understand the 
empirical correlation between electronic structure pa- 
rameters - the range of the in-plane hopping - and 
superconductivity, pointed out by Pavarini et al. Q? 
In particular, can we get some insight into the all im- 
portant question of what material parameters control 
the optimal SC transition temperature T™^^1 
Model: The minimal model that allows for an under- 
standing of the material dependencies of the cuprate 
pheonomenology is the t-J model with extended 
hopping: U = E^,j> (cLcjt + h. c.)P + 
J ' '^j ^ n-iUj/A) where Cia- is the electron 

operator at site i with spin a, riia = cj^Cia is the den- 
sity, with Hi = J^cr'^i'^^ = ^cj^^a/scip is the spin 
at site i (cr's are the Pauli matrices), and J is the an- 
tiferromagnetic exchange between nearest neighbors 
(i, j). The projection operator P = Yii ~ '^iT'^u) 
implements the "no double occupancy" constraint. 

The bare dispersion has the form e(k) = 
— 2i(cos kx + cos ky) + it' cos kx cos ky — 2t"(cos 2kx + 
cos2fcy) where t, t' and t" are the nearest, second 
and third-neighbor hoppings respectively. The im- 
portance of t' and t" is suggested both by ARPES 
experiments 0, S| and electronic structure calcula- 
tions 0. With the sign convention above, t,t' and 
t" are all positive for the hole doped case. To model 
the electron-doped case, we make a standard particle- 



hole transformation q = (— l)*Ci and cj = (—1) 
Ti. Thus for the electron-doped case we again obtain 
n with i = t, i' = -t', t" = -t". 
Variational Wave function: We choose a varia- 
tional ground state wave function for an A'^-particle 
system that includes both AF and SC order: 



(1) 



The form of in the unprojected wave function is 
motivated by a saddle point analysis of Ti. For a 
nonzero Neel amplitude m^, wc get two spin density 
wave (SDW) bands (a = 1,2) £;(i,2)(k) ^ [C(k) + 
C(k + Q)]/2 ± (l/2)v/(e(k) - e(k + Q))2 + WJ^mj, 
where Q ~ (tt, tt) and ^(k) = —2t{coskx + cos ky) + 
4i(,3^j, cos kx cos ky — 21'^^^ (cos 2kx +cos 2ky ) — /ivar , with 
k in the reduced Brillouin zone. The c?-wave pairing 
field A(k) = JA{coskx — cos ky) gives rise to SC in 
the two SDW bands, with BCS coherence factors 



A(k) 



E^{k) + Vi?^(k)+A2(k) 



(2) 



The internal pair wave function ip{ri — Vj) is given 
by the Fourier transforming Wka/ttka, and summing 
over the two SDW bands a. 

The five variational parameters in |\l/o) are the 
Neel amplitude rriN, the d-wave gap A, the (Hartree- 
shiftcd) chemical potential /.tvar, and t'^^,. and t'^^j. 
renormalized by Fock shifts. Their optimum values 
are determined by minimizing the ground state en- 
ergy E = ('I'o|W|*o)/(*o|*o) calculated using vari- 
ational Monte Carlo (VMC), which exactly imple- 
ments the projection P. We have developed a fast 
conjugate gradient algorithm (details will be pub- 
lished elsewhere) which evaluates derivatives of the 
energy and efficiently determines the optimized vari- 
ational parameters. 
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FIG. 1: (color online) (a) Phase Diagram of extended t-J model with t" = and t' indicated in the inset, showing 
Neel magnetization M and the SC order parameter $ as a function of both hole- and electron-doping. (b),(c) Shapes 
of the bare Fermi surface at half filling. 



Results: The T = phase diagram is deter- 
niincd by computing the SC and AF order param 
eters for the optimized ground state as a func- 



tion of doping [lO[. The SC order parameter [ll| 



lim 



cFa.ai'i' — f') is obtained from 
the long range behavior of the correlation func- 



|r-r' 



r') 



where, B. 



tion Fa^f}{r 

^{cr\'^r+ai~'^ri'^r+ai) creates a singlet on the bond 
(r,r -I- d), a = x,y. The AF order parameter 
M = (2/iV)(E.,^5f - EiesS!) is the difference 
of the magnetization on the A and B sublattices. 

In Fig. [1] we see the following phases: an AF Mott 
insulator at half-filling; coexistence of SC and AF 
for small (electron/hole) doping; a d-wave SC at 
higher doping; and a Fermi liquid for sufficiently large 
doping. The phase diagram shows marked electron- 
hole asymmetry for t' ^ 0. As the next-neighbor 
hopping \t'\ increases, SC is enhanced on the hole 
doped side, while AF is stabilized on the electron 
doped side. These results are consistent with earlier 
VMC 0, [H, m, O, M, m and dynami- 



2fl, 21^ 



cal mean field studies 

Our new findings are that, on the hole doped side, 
^'1 does not affect AF, and in particular the critical 
doping xaf, beyond which AF vanishes, is insensitive 
to the value of \t'\. On the electron doped side, SC is 
slightly weakened and the peak value of the SC order 
parameter ^max falls with increasing \t'\. 

Upon adding a second neighbor hopping t" wc find 
the following general trends: (a) On the hole doped 
side, an increase in t" leads to an increase of su- 



perconducting correlations. Interestingly, AF is rela- 
tively unaffected, and xaf is quite insensitive to t". 
(b) On the other hand, on the electron doped side, AF 
is enhanced with increasing \t"\, and superconductiv- 
ity is essentially unaffected. We have also performed 
calculations using a simple renormalized mean-field 
theory and the qualitative phase diagram is 

in agreement with the VMC results. 

At very low doping, our results differ from experi- 
ments due to the neglect of long range Coulomb and 
disorder effects. Once the local supcrfluid density, 
becomes sufficiently small near half-filling due to pro- 
jection Uj, long-wavelength, quantum phase fluctua- 
tions, neglected in our approach, drive SC to zero at a 
finite doping [i^l . We also do not consider the effects 
of other competing orders (stripes or charge order- 
ing). For large doping, where the pairing A — > 0, 
finite size effects become large because of a gro wing 
correlation length. This leads to an overestimate 24 1 
of the range over which SC survives, but does not 
qualitatively affect our conclusions. 

"Fermi surface" topology: Is there a simple way 
to understand our results? To begin with, we empha- 
size that the dependence of the phase diagram on the 
bare dispersion is not controlled by the van Hove sin- 
gularity (vHS) in the bare density of states. Even in 
the particle- hole symmetric case t' = t" = 0, where 
the vHS is precisely at the center of the band, the 
calculated SC order parameter $(x) does not peak 
at X = but rather at an 'optimal' doping away from 
zero. This 'optimal' doping is determined by the in- 
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FIG. 2: (color online) Left: A contour plot of as a 
function of the hopping amplitude t'/t and t" /t. Right: 
Monotonic relationship between 77 and the range param- 
eter 01 t'/t (for the special case of t" = t' /2). 

tcrplay between the growth of the pairing amplitude 
A with underdoping and the suppression of phase co- 
herence by strong correlations as a; ^ 0. Irrespective 
of the values of t' and t" , the renormalized bandwidth 
is of order {xt + J), while the scale of the pairing is 
also of order J. Thus we are not in a weak coupling 
BCS regime where all the action is in the immedi- 
ate vicinity of the chemical potential. In fact, the 
entire band participates in pairing, and the proxim- 
ity of the vHS to the chemical potential is not the 
dominant factor in determining the phase diagram. 

We next define a quantity 77 that characterizes the 
topology of the "Fermi surface" at half filling, defined 
by the bare tight-binding parameters. We show that 
this single parameter rj controls the dependence of the 
phase diagram on t' and t" . In particular, key fea- 
tures such as the critical doping xaf, the maximum 
value of the SC order parameter ^max, and optimal 
doping, are all determined by 77. 

For the cases of our interest, the "Fermi surface" 
at half filling can be described as a curve in the first 
Brillouin zone kpiO) (in polar coordinates) with 9 
measured from the fcr axis. We define 



77 {t'/t, t"/t) = 2 [kF{T:/A)/kF{emiu)Y 



(3) 



where ^min is the minimum angle at which a Fermi 
crossing exists. In the electron-doped case (FiglTJb)), 
Smin = and 7; > 1 corresponding to a convex "Fermi 
surface". In the hole-doped case (Fig[T]Jc)), 9min cor- 
responds to the Fermi crossing on the zone boundary 
leading to a concave "Fermi surface" with 77 < 1 . For 
t' = t" = Q, the particle-hole symmetric case, ?/ = 1. 
In Fig. [2] we show the dependence of 77 on t' and t" . 
The right panel focusses on the i'-dcpendcncc of rj on 
the special line t" = t' /2 which corresponds to the 
parameters obtained by Pavarini et al. [6| from elec- 
tronic structure calculations of single layer cuprates. 

In Fig. [3] we show the 77-dependence of the critical 
doping xaf for the vanishing of AF order. On the 
electron doped side, xaf increases approximately lin- 
early with 77, while for hole doping xaf is essentially 
independent of 77. Turning now to superconductiv- 
ity, we sec from Fig. [3] that the SC order parameter 
^max at optimality increases roughly linearly with 
decreasing rj on the hole doped side. On the electron 
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FIG. 3: (color online) Top: The size of AF phase xaf 
as a function of 77. Note that systems with very different 
bare ts but with the same 77 have similar xaf- This is 
also true for $max and optimal doping Xopt- Bottom: SG 
order at optimal doping $max as a function of 77. 
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FIG. 4: (color online) Left: Comparison of kinetic en- 
ergy per site obtained from variational Monte Carlo and 
Gutzwiller projected bare kinetic energy (BKE); {g = 
22;/(l -I- 2;)) Right: Bare kinetic energy as a function of 
doping for systems with different 77. 

doped side, there is a slight linear decrease of ^max 
as a function of increasing rj. Thus a more concave 
bare "Fermi surface" leads to a more stable SC state. 
We have also studied the dependence of the "optimal 
doping" Xopt (doping at which $ attains ^max) and 
found a similar correlation with 77. 

We have thus demonstrated that three character- 
istics of the phase diagram, xaf, ^max and Xopt, are 
determined by a single parameter rj, rather than by 
the details of the bare dispersion. Two systems with 
a given t and J but with two very different t', t" have 
the same phase diagram provided they correspond to 
the same value of 77. This point is clearly illustrated 
by the specially marked points in Fig. [3] 

To understand how the single parameter rj con- 
trols the entire phase diagram, we consider the com- 
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petition between the kinetic and exchange energies. 
Upon doping the Mott insulator, the holes (or elec- 
trons) attempt to gain kinetic energy (KE). Near half 
filling, the hopping t between different sublattices dis- 
rupts the anti-fcrromagnctic order and increases the 
exchange energy. On the other hand, the hoppings 
t' and t" between the same sublatticc do not disturb 
AF order and entail no exchange energy penalty. In- 
sight into how the system gains the most KE, while 
keeping the exchange energy increase to a minimum, 
may be obtained from studying the behavior of the 
KE as a function of doping x. 

In Fig. BJa), we show that the variational Monte 
Carlo (projected) KE for a convex (77 > 1), a con- 
cave (77 < 1), and a half-filled diamond (77 = 1) Fermi 
surface. In each case the projected KE is closely re- 
produced by the Gutzwillcr approximation [l9| result 
g Kbare whcrc Kbarc IS thc bare KE in the unprojected 
state, and g = 2x / {1 + x) the rcnormalization factor 
that takes into account projection. 

The renormalization factor g is independent of the 
77-parameter, and thus we focus on the bare KE in 
Fig. Ub) to understand the 77-dependence of the com- 
petition between KE and superexchange. For elec- 
tron doping with 77 > 1, the bare KE is a non- 
monotonic function of doping with a minimum at a 
finite Xm- Thus up to a doping of Xm, one can gain 
KE due to t' and t" without sacrificing exchange en- 
ergy, thereby stabilizing the AF state. The doping 
Xm increases with increasing rj (see Fig. |4l[b)) which 
underlies the r^-dependence of xaf on the electron 
doped side. For the p-h symmetric {t' = t" = 0) 
case, 77 = 1 and the minimum KE is at Xm = 0. For 
hole doping with rj < 1, the bare KE increases mono- 
tonically with doping. The exchange energy is best 
satisfied by means of SC with resonating singlet pairs, 
while providing for the necessary KE gain, respecting 
the no double occupancy constraint. Thus a system 
with a larger bare KE favors a more stable SC state. 

As noted earlier, Pavarini et al. 0] suggested an 
empirical correlation between the range parameter 
related to t' /t (for the special case of t" = \t'\/2) and 
j-max^ the maximum Tc within a given cuprate fam- 
ily. We cannot, of course, obtain from our ground 
state calculation, but the magnitude of the SC order 



parameter $ can be taken [ll|, |17| as a measure of 
Tc- The strong linear correlation between ^"^"^ and 
r] on the hole-doped side seen in Fig. [3] is equivalent 
to a similar correlation between T™'^^ and the range 
parameter. (Note the relationship between rj and the 
range parameter in Fig. [2]). On the electron-doped 
side, however, we do not predict a strong dependence 
of T^^^ on thc 7] parameter. 

In conclusion, our results provide a unified micro- 
scopic picture for the competition between AF and 
SC for electron and hole doped cuprates, while also 
providing insights into a key material parameter 77 



that controls many aspects of the phase diagram, in- 
cluding T™''^. In particular, this work suggests a 
route to increase T™°^ by creating systems with a 
small 77, i. e., with a highly concave bare "Fermi sur- 
face" at half filing. Even if such band structure en- 
gineering may not be easy in solid state materials, 
it may be possible in optical lattice realizations of 
strongly correlated Fermions. 
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